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ABSTRACT 

The results from deep long-slit spectroscopy of nine edge-on spiral galaxies with known extraplanar line emis- 
sion are reported. Emission from Ha, [N II]AA6548, 6583, and [S II]AA6716, 6731 is detected out to heights of a 
few kpc in all of these galaxies. Several other fainter diagnostic lines such as [O I] A6300, [O III] AA4959, 5007, 
and He I A5876 are also detected over a smaller scale. The relative strengths, centroids and widths of the various 
emission lines provide constraints on the electron density, temperature, reddening, source(s) of ionization, and 
kinematics of the extraplanar gas. In all but one galaxy, photoionization by massive OB stars alone has difficulties 
explaining all of the line ratios in the extraplanar gas. Hybrid models that combine photoionization by OB stars 
and another source of ionization such as photoionization by turbulent mixing layers or shocks provide a better 
fit to the data. The (upper limits on the) velocity gradients measured in these galaxies are consistent with the 
predictions of the galactic fountain model to within the accuracy of the measurements. 

Subject headings: diffuse radiation - galaxies: halos - galaxies: ISM - galaxies: spiral - galaxies: structure 
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1. INTRODUCTION 

Emission-line diagnostics have been used successfully to de- 
termine the hardness of the ionizing spectrum in Galactic and 
extragalactic H II regions (e.g., Stasinska 1982; Evans & Dopita 
1985; McCall, Rybski, & Shields 1985; Dopita et al. 2000 and 
references therein) and in the nuclei of galaxies (e.g., Baldwin, 
Phillips, & Terlevich 1981; Veilleux & Osterbrock 1987; Os- 
terbrock, Tran, & Veilleux 1992; Veilleux 2002 and references 
therein), but only over the last decade has it been possible to 
measure the emission line ratios in the faint, extraplanar dif- 
fuse ionized gas (eDIG) of external galaxies (e.g., reviews by 
Dettmar 1992 and Dahlem 1997). Observations of the diffuse 
ionized gas in our own Galaxy show line ratios which are dif- 
ficult to explain with pure stellar photoionization models with- 
out extra heating (e.g., Reynolds 1985a, 1985b; Reynolds & 
Tufte 1995; Mathis 2000). A similar situation appears to apply 
to external galaxies. The [N II] A6583/Ha and [S II] A6716, 
6731/Ha line ratios measured in a few galaxies generaly be- 
come stronger with increasing heights, often reaching values 
considerably higher than typical values observed in H II re- 
gions (e.g.. Rand, Kulkarni, & Hester 1990; Keppel et al. 1991; 
Dettmar & Schultz 1992; Veilleux, Cecil, & Bland-Hawthorn 
1995; Ferguson, Wyse, & Freeman 1996; Golla, Dettmar, & 
Domgorgen, 1996; Domgorgen & Dettmar 1997; Rand 1998; 
Otte & Dettmai- 1999; Tullman & Dettmar 2000; Tullman et al. 
2000; Miller & Veilleux 2003a, hereafter Paper I). The verti- 
cal [N II]/Ha and [S II]/Ha gradients detected in these galaxies 
may be due to hardening of the OB-star radiation as it passes 
through the dusty and neutral medium of the galaxy, or to the 
existence of other sources of heating or ionization which is be- 
coming increasingly important above the galactic plane. Pos- 
sible sources of extra ionization and heating include shocks, 
photoionization by cooling hot gas, "turbulent mixing layers" 
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(TML; Slavin, ShuU, & Begelman 1993) or supernova remnants 
(Slavin, McKee, & Hollenbach 2000), cosmic ray heating (e.g., 
Lerche & SchHckeiser 1982; Hartquist & Morfill 1986; Parker 
1992), and magnetic reconnection (e.g., Raymond 1992). 

The measurements of additional line ratios can shed some 
light on the importance of secondary ionization sources. One 
particularly important line ratio is [O III] A5007/Ha, a good 
indicator of high energy processes. However, [O III] is chal- 
lengingly faint and has therefore been measured in only a few 
galaxies (e.g.. Rand 1998; Tullman & Dettmar 2000; Tullman 
et al. 2000; Collins & Rand 2001). This is also the case for He 
I A5876/Hq!, a sensitive indicator of the hardness of the ioniz- 
ing radiation. Interestingly, the value of He I/Ha in NGC 891, 
a galaxy which in many ways is very similar to our own, ap- 
pears significantly larger than the Galactic value (0.034 versus 
0.012 ± 0.006; Reynolds & Tufte 1995; Rand 1997), while the 
value observed in NGC 3044 is even larger 0.07; Tullman 
& Dettmar 2000). [O II] A3727/H/? has recently been shown to 
be a useful diagnostic of extra heating in the diffuse ionized gas 
(Mathis 2000; Otte et al. 2001; Otte, Gallagher, & Reynolds 
2002), but li/3 emission is generally very faint outside of H II 
regions and [O II] A3727/Ha is highly sensitive to reddening 
corrections and flux calibration errors. 

There is a need to expand the set of high-quality spectro- 
scopic observations of the eDIG to a larger number of galaxies. 
This paper describes an attempt to remedy this situation. The 
results from a spectroscopic survey of nine edge-on galaxies 
with known eDIG are reported. Due to scheduling constraints, 
the imaging observations reported in Paper I were not reduced 
and analyzed in time for our scheduled spectroscopic observa- 
tions, so the spectroscopic sample was selected independently 
of the imaging sample. The only exception is NGC 2820, where 
the Ha image obtained with the TTF (see Paper I) was used to 
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determine the optimum slit position [note that NGC 4013 is also 
in the imaging sample of Paper I, but the position of the slit for 
this object is based on earlier observations by Rand ( 1 996; here- 
after R96)]. The nine galaxies in the spectroscopic sample were 
selected based on the published reports of extraplanar emission 
by Pildis, Bregman, & Schombert (1994b, hereafter PBS) and 
R96. The methods used to acquire and reduce these data are 
discussed in §2. Great care is taken to reach a limiting sur- 
face brightness of order a few 10"^^ erg s~' cm~^ arcsec"^. The 
results from the spectroscopic analysis are given in §3. New 
constraints on the physical conditions in the eDIG (e.g., tem- 
perature, density, reddening, ionization level, kinematics) are 
derived using the relative strengths and positions of the stronger 
emission lines that lie within 4550 - 7300 A. In §4, the line 
ratios derived from the long-slit spectra are compared with val- 
ues measured in other galaxies, as well as with predictions from 
photoionization models (Sokolowski 1994; Bland-Hawthorn et 
al. 1997), turbulent mixing layer models (Slavin et al. 1993), 
and shock models (ShuU & McKee 1979; Dopita & Suther- 
land 1995). This analysis allows us to determine whether a sec- 
ondary source of ionization in addition to photoionization by 
hot stars is needed in the eDIG. The main results are summa- 
rized in §5. 

2. DATA ACQUISITION AND REDUCTION 

AH of the data were taken at the KPNO 2.1-m telescope on 
January 22 - 26, 1998. A compromise had to be made between 
broad wavelength coverage and good spectral resolution. Grat- 
ing #26new and filter GG-420 were used with the F3KA CCD 
to provide a dispersion of 1.256 A pixel"' and a useful spectral 
coverage of ~ 2760 A between ~ 4550 and 7300 A, after ac- 
counting for the known bad columns on the redward side of the 
F3KA CCD. The slit width was set at 1"5, yielding a spectral 
resolution of 3.7 A. This configuration allows us to resolve 
the Ha + [N H] AA6548, 6583 complex, the [S H] AA6717, 673 1 
doublet, and the He I A5876 line from the nearby Na ID sky 
lines, but does not cover the important [O II] A3727 and [O III] 
A4363 diagnostic lines. The length of the slit in this mode is 
5 '.2, extending well beyond the extent of the galaxy disks in our 
sample. The CCD was binned in the spatial direction by a fac- 
tor of 2 in order to increase the signal per pixel from the diffuse 
emission, resulting in a spatial scale of 1 "56 pixel"' . To reach 
our goal of achieving a flux limit on the order of a few x 10"'^ 
erg s~' cm"^ arcsec"^, each galaxy was observed for about 5 
hours. Details on the observations are listed in Table 2. The 
slit was centered on the disk of the galaxy and in most cases 
positioned so that it lay perpendicular to the disk, although in a 
few cases it was tilted slightly to optimize the coverage of the 
extraplanar emission. The position of the sUt is mentioned in 
Table 2 and shown in Figure 1 for each galaxy in the sample. 
All galaxies were observe through an airmass of less than ^ 1 .5 
to avoid any significant differential atmospheric refraction. 

Bias frames were taken each night and a composite bias was 
made by combining the individual frames. The bias level was 
found not to be constant across the CCD, showing a gradient 
along the dispersion axis. A one-dimensional fit was applied 
along the dispersion axis to create the bias frame, making sure 
not to introduce additional noise to the data when subtracting 
off the bias. The spectra were then trimmed and corrected for 
bias and overscan using the CCDPROC package in IRAF. Both 
dome and internal (quartz) flats were obtained during the night 
with the purpose to use them for the flatfield and illumination 



corrections. However, better results were obtained when us- 
ing the domeflats to flatfield the data, and the data themselves 
to correct for the illumination along the slit. For this procedure, 
each spectrum was binned along the dispersion axis so that each 
wavelength bin was well represented and contained sufficient 
counts. The illumination variations across the slit were accu- 
rately modeled using the IRAF ILLUMINATION routine to fit 
a spline function to the background while avoiding the emission 
due to the galaxy. 

Next the data were corrected for distortion and wavelength 
calibrated using HeNeAr spectra. Observations of the HeNeAr 
lamp were obtained before and/or after each galaxy observation 
so that accurate corrections could be obtained. A fit to the sky 
background was calculated and subtracted from each frame us- 
ing sample rows on either side of the galaxy, far enough from it 
to avoid the extraplanar emission. The two-dimensional spec- 
tra were then flux calibrated using HZ 44 as a calibrator and 
combined together. 

After combining the different observations for each galaxy, 
the redshifted emission lines were identified based on the sys- 
temic optical velocity of each galaxy (as listed in de Vau- 
couleurs et al. 1991). The spatial dimension was further binned 
by two in order to increase the signal per binned pixel in the 
weaker fines, and one-dimensional spectra were extracted at 
different heights from the disk plane. The spatial dimension 
was binned even further at large |z| to help detect very faint ex- 
tended line emission. The spectra for each galaxy are shown 
in Figure 2 as a function of position along the slit. In the case 
where a spectral line is not detected, a representative spectrum 
from near the disk plane is presented. The line fluxes were 
determined using the SPLOT routine in IRAF. Single gaussian 
profiles were fitted to each fine in order to calculate the fine flux 
and width and the level of the underlying continuum emission. 

3. RESULTS 

3.1. Line Ratios 

The vertical profiles of up to nine diagnostic line ratios are 
shown in Figure 3 for each galaxy. In each figure, the vertical 
profile of the Ha emission is shown as the heavy solid fine and 
the continuum emission is shown as the dotted line. Both have 
been scaled arbitrarily for display purposes. To prevent the fig- 
ures from becoming overcrowded, the error bars in all of the 
figures represent Icr uncertainties as determined from SPLOT. 
In most galaxies. Ha, [N II] A6583, and [S II] AA6716, 6731 
have been detected along the sHt out to 1 - 2 kpc from the center 
of the disk (in most of our galaxies, the slit is perpendicular to 
the disk and this represents the actual vertical extent of the gas, 
but in a few cases this distance differs slightly from the vertical 
height; see Fig. 1). Fainter lines such as H/?, [O III] A4959, 
5007, and He I A5876 are also detected over a smaller scale in 
a number of galaxies. This section of the paper discusses the 
overall trends found in the sample. For a more detailed discus- 
sion of each object, the reader should refer to the Appendix. 

3.1.1. [N II]/Ha, [S II]/Ha, and [N II]/[S II] 

The average midplane values for [N II] A6583/Ha and [S II] 
A6716/Ha are 0.40 ± 0.20 and 0.24 ± 0.14, respectively. Seven 
galaxies show a general increase in [N II]/Ha and [S II]/Ha 
with increasing height. The most dramatic gradients occur 
wifliin NGC 4013 and NGC 4217, where the [N II]/Ha ratios 
are H II region-Uke in the disk (~ 0.3 - 0.4) but reach val- 
ues of nearly 2.0 at |z| 1 - 2 kpc above the disk. Some of 
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the galaxies in the sample have [S II] A6716/[N II] A6583 ra- 
tios which change with height; some (such as NGC 2820 and 
UGC 4278) show a steady increase with increasing |z|. Others 
(such as NGC 4302) show a general decrease with increasing 
height. Most interesting is NGC 4217, which shows an increase 
in [S II]/[N II] from 0.5 to 0.7 up to |z| ~1 kpc, and then the ra- 
tio falls to 0.3 at higher |z|. Other galaxies present [S II]/[N II] 
ratios which are consistent with being constant within the un- 
certainties of the measurements. 

For comparison, Rand (1998) found in NGC 891 that the 
[N II]/Ha line ratio rises from 0.35 in the plane of this galaxy, to 
a value greater than unity at |z| ^ 2 - 3 kpc. A similar trend was 
found with the [S IIJ/Ho; Une ratio, such that the [S II]/[N II] 
ratio remained almost constant with a value ~ 0.6. Observa- 
tions of a few other galaxies (Collins & Rand 2001; Tiillmann 
& Dettmar 2000) find similar trends. The ratio of coUisionally 
excited lines (Uke [N II] and [S II]) to recombination lines (like 
Ha) depends on the ratio of heating to recombination. Since 
S° and N" have similar ionization potentials (10.4 eV and 14.5 
eV, respectively), any changes in the |S 11|/|N II] ratio suggest 
changes in the local ionization condition (S^ ionizes at a slightly 
lower ionization energy than N"""), or in the metallicity (since ni- 
trogen is a secondary product of nucleosynthesis while sulphur 
is a primary product). These issues are discussed further in §4.2 
and §4.3 below. 

3.1.2. He I/Ha 

He I A5876 is detected in the disk of only 4 of the 9 sam- 
ple galaxies (the presence of the Na ID absorption lines in at 
least 4 of the galaxies makes detection of this line impossible 
- see Fig. 2 for details). Emission from H 11 regions in the 
disks of the galaxies is almost certainly contaminating some of 
these measurements. The midplane value of He 1/Ha for these 
galaxies ranges from 0.018 (NGC 4013) to 0.052 (NGC 4302), 
bracketing the value measured in Orion (0.042; Kaler 1976) but 
larger than the value measured in the diffuse ionized gas near 
the midplane of our Galaxy (0.012 ± 0.006; Reynolds & Tufte 
1995). The relative ionization fractions of helium and hydrogen 
can be determined by the equation: 
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(e.g., Brockelhurst 1971; Martin 1988; Reynolds & Tufte 1995; 
Rand 1997), where E is the emissivity in cm~^ s~^ , % is the frac- 
tion of helium or hydrogen that is singly ionized, He/H is the 
abundance of helium with respect to hydrogen by number, and 
T is the gas temperature. Using the He/H abundance listed in 
Boesgaard & Steigman (1995) and J = 8000 K (e.g., Reynolds 
1992), xhJxh is found to range from 0.38 to 1 .09. If hydrogen 
is assumed to be mostly ionized (a reasonable assumption given 
the strength of [O III] in the eDlG of several galaxies; see §4.2), 
then the helium is about 40% ionized in NGC 4013, and is al- 
most fully ionized in NGC 4302. Using Table 1 of Rand (1997), 
QhcIQh, the ratio of He-ionizing {hv > 24.6 eV) to H-ionizing 
{hv > 13.6 eV) photons, ranges from ^ 0.040 to 0.115. These 
results imply an effective temperature of the radiation field, J*, 
which ranges from about 36,500 to 38,500 K, and an upper limit 
to the stellar mass function, M„, which ranges from ~ 42 to 54 
Mq. 

Extraplanar He I emission is detected unambiguously in only 
one galaxy of our sample, NGC 2820. This object has a mid- 
plane He I/Ha value of 0.046, and a range of values from 0.033 



to 0.056 (the latter at a height of |z| = 1 kpc). Correspond- 
ing values of [N II]/Ha are ^ 0.25 at the midplane, and range 
from 0.16 to ~ 0.3 at |z| = 1 kpc. Within the uncertainties, 
these values are consistent with the predictions from the O-star 
photoionization models of Domgorgen & Mathis (1994). Rand 
(1997) had difficulties reconciling the values of He I/Ha with 
those of [N II]/Ha in NGC 891, but the extraplanar [N II]/Ha 
line ratios in this galaxy are much higher (near 1.4). Extrapla- 
nar He I was not detected in any of the galaxies of our sample 
with large extraplanar [N II]/Ha ratios. 

3.1.3. Reddening from Ho/HP 

The Ha/H/3 line ratio is shown in the bottom left panel of 
Figure 3 for each galaxy (when detected). No correction was 
made for possible underlying stellar absorption features; these 
features were never evident in the data. This line ratio provides 
an indication of the amount of reddening affecting the spectra. 
Determining the amount of reddening is not straightforward, 
however. Typical Ha/H/3 line ratios for H II regions is ^ 2.85, 
but if shocks are present, this ratio could rise to ~ 3.10 (e.g., 
Shull & McKee 1979). The effects of extinction and reddening 
also depend on the distribution of the dust with respect to the 
source of emission (e.g., uniform screen in front of the source 
of emission versus uniform dust distribution mixed with the line 
emitting gas). Rather than to try to correct for these complex 
effects, we instead list in Table 3 the impact on the line ratios 
of a foreground screen of dust with Ay = 1. Except for Ha/H/3 
and [O III]/Ha, the line ratios shown in Figure 3 are not at all 
sensitive to reddening. 

In general the amount of reddening is found to be larger in 
the disk than in the eDIG, as one would expect if the dust is dis- 
tributed near the plan of the galaxy disk. Out of the five galaxies 
in which extraplanar H/3 was detected, three clearly show this 
trend. In the other two galaxies, the Ha/H/3 line ratio appears 
to be relatively constant with height (NGC 2820) or does not 
show any obvious monotonic trend with height (NGC 4013). 

3.1.4. Density from [S II] X6716/[S II X6731 

The [S II] A6716/[S II] A6731 ratio is shown in the bottom 

middle frame of Figure 3 for each galaxy (when detected). The 
ratio of the intensities of these lines yield information on the 
average electron density of the gas. The low density limit of 
this line ratio is ~ 1.4, and so it is possible to make quantita- 
tive statements regarding the density of the gas only in regions 
where this line ratio is less than ^ 1.4. 

For most of the galaxies, the [S II] A6716/[S II] A6731 line 
ratio is consistent with the low density Umit, therefore suggest- 
ing an electron density of at most a few tens of cm"^. This 
is typical of the electron density that has been reported for 
other galaxies (e.g.. Rand 1998; Collins & Rand 2001). There 
are three exceptions: NGC 3628, NGC 4217, and NGC 4302. 
The results on NGC 3628 are discussed in the Appendix. In 
NGC 4217, the [S H] A6716/[S II] A6731 line ratio has a value 
of ~ 1.2 near the disk of the galaxy, and drops slightly to an 
average value of 0.9 at a height of about 1 kpc, before climbing 
back up over 1.4. Assuming a constant temperature of ~ 10"* 
K for the moment, this suggests that the electron density in the 
disk of the galaxy is ~ 200 cm~^ and increases to a value of 
900 cm"^ at heights ~ 1 kpc, before decreasing below the low 
density Hmit. In NGC 4302, the value of the [S II] A6716/[S II] 
A6731 Une ratio near the disk of the galaxy is ~ 1.0, and de- 
creases to values near 0.8 and 0.6 at heights of 0.5 kpc. The 
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corresponding electron densities are ^ 500 cm~^ in the disk of 
the galaxy and 1000 to 2000 cm~^ at larger heights. It should 
be noted, however, that the electron density as measured by the 
[S II] ratio scales as T^/^ (e.g., Osterbrock 1989), and therefore 
an increase in temperature would be interpreted as higher elec- 
tron density if constant temperature were assumed. The vertical 
temperature profiles in the eDIG of these galaxies are discussed 
next. 

3.1.5. Temperature from [N IIJ/Ha 

Two of the best temperature gauges for the extraplanar gas 
are the [N II] A6583/[N II] A5755 and [O III]A5007/[O III] 
A4363 line ratios. [N II] A5755 was recently detected in the 
diffuse ionized gas of our Galaxy and indicates elevated temper- 
atures relative to those of H II regions (Reynolds et al. 2001). 
Unfortunately the [N II] A5755 line was not detected in any of 
the sample galaxies and [O III] A4363 falls outside the wave- 
length range of our observations. We therefore have no choice 
but to use other Une ratios for this analysis. Recent studies have 
suggested the use of the [N II]/Ha line ratio as a temperature 
diagnostic (Haffner, Reynolds, & Tufte 1999), keeping in mind 
that the derived temperature is an average value over the line- 
of-sight column density. Given that hydrogen and nitrogen have 
similar first ionization potentials and a weak charge-exchange 
reaction, and assuming that hardly any N is doubly ionized, one 
has N'^/N*' « H^/H". Under this assumption and using a Galac- 
tic gas-phase abundance of (N/H) = 7.5 x 10"^ at all heights 
(Meyer, Cardelli, & Sofia 1997), the relationship between the 
[N II]/Ha Une ratio and electron temperature is given by 

^^ = I2.2rr2«e-2-i8/^^ (2) 
Ha 

where T4 is the electron temperature in units of lO'' K (ColUns 
& Rand 2001). The detection of significant [O III] emission in 
the eDIG of four galaxies of our sample ([O III]/Hq; > 0.25; 
UGC 4278, NGC 2820, NGC 4302 and NGC 5777) is incon- 
sistent with the [N II]/Ha temperature model (eqn. 2) because 
this model assumes that is not present in the eDIG, yet 
the ionization potential of (29.6 eV) is close to that of O""" 
(35.1 eV). Temperature profiles for the remaining five galax- 
ies in the sample have been calculated using eqn. (2) and the 
measured [N II]/Ha line ratios; the results are shown in the bot- 
tom, left panel of Figure 4 (the other panels will be discussed 
in §4.2). Three of these galaxies (NGC 3628, NGC 4013 and 
NGC 4217) show an increase in temperature with vertical dis- 
tance from the disk. The change in temperature with height 
varies from galaxy to galaxy, with some objects showing no 
change with height within the uncertainties (e.g., UGC 2092 
and UGC 3326), while others show a dramatic increase (e.g., 
NGC 4217: from 6000 K to almost 10,000 K). Haffner et al. 
(1999) detect an increase in electron temperature for our Galaxy 
from 7000 K at |z| = 0.75 kpc to over 10,000 K at |z| = 1 .75 kpc. 
Collins & Rand (2001) find similar temperature gradients in 
their studies of edge-on galaxies. 

3.2. Kinematics 

The kinematics of the extraplanar gas can also provide in- 
sights into its nature and origin. To our knowledge, only two 
galaxies (NGC 891 and NGC 5775) have so far been exam- 
ined with sufficient care to address this important issue. In 
both cases, the velocity of the extraplanar material is seen to 
approach the systemic velocity of the galaxy (Pildis, Bregman, 



& Schommer 1994a; Rand 1997, 2000; Tullman et al. 2000). 
This vertical velocity gradient can be explained by a combina- 
tion of radial movement of the gas (following the pressure gra- 
dient of the halo) and declining rotation speed (conserving an- 
gular momentum), as predicted by the galactic fountain model 
(e.g., Bregman 1980; Houck & Bregman 1990) 

Although our spectroscopic setup was chosen to optimize 
wavelength coverage at the expense of spectral resolution to al- 
low us to carry out a detailed line ratio analysis of our objects, 
an attempt is made to constrain the kinematics of the eDIG in 
the galaxies of our sample. Figure 5 displays the velocity pro- 
files for Ha, [N H] A6583, and [S U] A6716, the strongest Unes 
in our galaxies. These velocity profiles are not expected to dif- 
fer significantly from each other, so comparisons between the 
three panels can serve to estimate the uncertainties of the ve- 
locity measurements. Galactic rotation causes a velocity offset 
from systemic in the cases where the sUt does not pass through 
the nucleus. Figure 5 shows that most galaxies do not show 
significant vertical gradients within the accuracy of the mea- 
surements (~ 30 - 50 km s~' depending on the galaxy). How- 
ever, there appears to be a few exceptions. A significant dip 
in velocity {5v « 100 km s~') is observed at z « -0.3 kpc in 
NGC 3628, coincident with the prominent dust lane in this ob- 
ject. Dust obscuration is severely limiting the use of the optical 
emission lines as kinematic probes in this region. Possibly sig- 
nificant gradients are also visible in NGC 2820 and NGC 4013 
(and perhaps in NGC 4302 but at a lower significance level). 
In NGC 2820, the velocities reach a maximum of ^ 150 km 
s"^ around z = +0.4 kpc and then show a monotonic decrease 
toward systemic velocity at large heights, reaching values of ~ 
80 km s~^ at z w kpc. A slightly less significant gradient of 
~ 50 km s~^ appears to be present on both sides of the disk of 
NGC 4013 out to \z\ ~ 1 kpc. The gradients observed in both of 
these galaxies can be explained if the rotational velocity of line- 
emitting material is lower above and below the galaxy disk, as 
expected in the galaxy fountain model. The amplitudes of the 
detected gradients are however larger than the predictions from 
the model of Bregman (1980). The balUstic model of Collins, 
Benjamin, & Rand (2002) has more success explaining the gra- 
dients in NGC 2820 and NGC 4013. The lack of obvious gra- 
dients in the other galaxies is not inconsistent with the predic- 
tions of the galactic fountain model of Bregman (1980), given 
the relatively large uncertainties in the measurements. 

The widths of the emission lines provides another constraint 
on the gas kinematics. The lower right panels of Figure 5 
show the line widths of Ha, selected because it is generally the 
strongest emission line in our galaxies. Typical values he be- 
tween ~ 100 and 200 km s~^ There is Uttle or no evidence for 
significant line width gradients in the majority of the galaxies in 
our sample. The only galaxies where gradients may be present 
are NGC 3628 (coincident with the dust lane), UGC 3326 (the 
line widths vary from ~ 275 km s~^ in the midplane to ~ 150 
km s~^ at heights near |z| = 1 kpc), NGC 4013 (a sHght positive 
gradient of ^ 50 km s"' over ± 0.5 kpc on both sides of the 
disk may be present), and NGC 2820 (the line widths in this 
object are constant over most of the slit, but show a significant 
decrease from 170 km s"' at |z| ~ 1 kpc to ~ 80 km s"' at |z| ~ 
1.7 kpc in the southern halo of the galaxy). Numerical broad- 
ening arising from Poisson noise, whereby noisier signals tend 
to pull in more emission in the wings than at the peak, may af- 
fect (increase) the line widths at large |z|. Positive vertical line 
width gradients may also be due to increasingly turbulent mo- 
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tions at large heights (e.g., turbulent mixing layer; §4.3.1), or 
due to the fact that extinction by the disk is less significant at 
large heights (§3.1.3) and therefore a longer column of material 
with a broader range of velocities is being sampled. Our data do 
not allow us to distinguish between these various possibilities. 

4. DISCUSSION 

Several models have been proposed to explain the line ra- 
tios detected in the extraplanar material of disk galaxies. In 
this section we discuss each of these models and compare their 
predictions with our data. 

4.1. Photoionizationby OB Stars 

OB stars are almost certainly contributing to the ionization 
of the eDIG. They are by far the main source of Lyman contin- 
uum flux produced in the disk (e.g., Reynolds 1984), but their 
position near the disk midplane makes them highly vulnera- 
ble to absorption by the ISM and dust. Photoionization models 
(e.g., Mathis 1986; Domgorgen & Mathis 1994; Sokolowski 
1994; Bland-Hawthorn et al. 1997; Mathis 2000) have had 
some success explaining the increase in the [N IIJ/Ha line ratio 
observed in the eDIG. This increase is attributed to a decrease 
with height of the ionization parameter ([/), a measure of the 
ratio of the ionizing photon number density (<!') to the electron 
density (n,,). Under the assumption of ionization equilibrium, 
$ oc at all heights. Therefore, U oc oc so the ion- 
ization parameter should fall off exponentially with height. As 
U decreases, lower ionization species are favored, leading to an 
increase in the [N II]/Hq: and [S II]/Hq: ratios and a decrease in 
the [O III]/Ha ratio (neglecting the effects of reddening). 

However, some problems arise with the photoionization 
models. First, pure stellar continua models have difficulty re- 
producing [N II]/Hq; and [S Ill/Ha ratios greater than unity, as 
detected in some of the galaxies of our sample and in other stud- 
ies (e.g., Dettmar & Schultz 1992; Veilleux et al. 1995; Rand 
1998; Collins & Rand 2001; Paper I). Photoionization mod- 
els that take into account the multi-phase nature of the ISM, 
the possible depletion of certain gas-phase abundance of metals 
onto dust grains, and the absorption and hardening of the stel- 
lar radiation field as it propagates through the dust and H I gas 
in the disk (Sokolowski 1994; Bland-Hawthorn et al. 1997) are 
more successful at producing elevated [N Ill/Ha ratios of ~ 1 .5. 
Unfortunately, these models have difficulties explaining the ob- 
served behavior of [O IJ and [O III] relative to Ha. Collins & 
Rand (2001) detect [O III] in three out of four of their galaxies, 
and observe in each one a general increase in the [O III]/Ha 
line ratio with increasing height. Out of the five galaxies in 
which we detect [O III], three of them (NGC 2820, NGC 4302, 
and NGC 5777) show the same positive trend with height. A 
sharp increase in the [O I]/Ha fine, which is difficult to explain 
using photoionization models, is also observed in two galaxies 
of our sample (NGC 2820 and NGC 4217; see also Collins & 
Rand 2001; Rand 1998). A secondary source of heating and/or 
ionization appears to be needed to explain these observations. 

4.2. Secondary Source of Nonionizing Heating 

We first explore the possibility of a nonionizing source of 
heating in the eDIG (e.g., photoelectric heating from dust 
grains, dissipation of interstellar turbulence). The apparent 
temperature gradients found in three galaxies in the sample 
(§3.1.5) bring support to this scenario but do not prove it. For 
this, one needs to also examine the behavior of the other line 



ratios. Using the assumptions mentioned in §3.1.5 when deriv- 
ing equation (2), Galactic gas-phase abundances of S, N, and 
O, and the fact that the ionization fractions of oxygen and hy- 
drogen are coupled through a charge exchange reaction, Collins 
and Rand (2001; also Haffner et al. 1999) derive the following 
equations for the line ratio intensities as a function of tempera- 
ture and ionization fraction: 
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These equations imphcitly assume that the metallicity of the 
eDIG is constant with height. Note that eqn. (4) predicts only 
a very weak temperature dependence for the [S II]/[N II] ratio. 
Using the temperature gradients derived in §3.1.5, the line ra- 
tios can be predicted for a number of ionization fractions and 
compared with the measured quantities; this is done in Figure 4. 
The variations of the [S II]/[N II] line ratios in UGC 3326, 
NGC 3628, NGC 4013, and NGC 4217 (the large uncertainties 
in the data of UGC 2092 prevent us from making any statement 
for this object) cannot be simply explained by a temperature 
gradient with height. The abundance of N"^ relative to that of 
S"^ varies with heights, perhaps an indication that an additional 
source of ionization is present in these objects. 

The positive gradients in the vertical [O 111]/Hq: profiles of 
NGC 2820, NGC 4302, and NGC 5777 sfi-ongly suggest the 
presence of an additional source of ionization that becomes 
more influential at lower densities. This possibility has also 
been suggested to explain the line ratios in our Galaxy (Haffner 
et al. 1999; Reynolds et al. 1999) and a few other external 
galaxies (Collins & Rand 2001). The next section considers 
this scenario in more detail. 

4.3. Secondary Source of Ionization 

A wide variety of secondary sources of ionization have been 
suggested to account for the line ratios in the eDIG (see refer- 
ences in §1). In this section we discuss the predictions from 
three of the possible scenarios and compare them to our data. 
The promising possibility that cooling supernova remnants also 
contribute to the ionization of the eDIG (Slavin, McKee, & Hol- 
lenbach 2000) is not discussed here because the predicted line 
ratios are not available in this case. 

4.3.1. Turbulent Mixing Layers 

Slavin et al. (1993) model the turbulent mixing layer (TML) 
created by shear flows between hot and cold gas. This layer 
is made up of intermediate-temperature gas which radiatively 
cools until the cooling rate is balanced by the energy flux into 
the gas layer. The properties of the TML depend on two main 
parameters, namely the temperature attained by the gas imme- 
diately after mixing and the enthalpy flux per particle into the 
layer (which itself depends on the velocity of the hot gas, V(). 
Their models span a temperature range of 5.0 < log T <5.5 and 
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a velocity range of 25 km s"' < v, < 100 km s~'. In their mod- 
els, Slavin et al. assume pressure and ionization equilibrium 
between the hot and cold gas phases and slow mixing within 
the turbulent layer, so as to distinguish between turbulent mix- 
ing and shocks. 

In general, the TML models are successful at predicting 
many of the observed line ratios, such as [N II]/Ha and 
[S II]/Ha. The TML models also predict strong [O 111] emission 
(due to the high level of ionization in the mixing layer region), 
an emission line which is difficult to explain with photoioniza- 
tion alone (§4.1). 

Figure 6 presents plots of the [S nj/Ho!, [O m]/Ha and 
[O I]/Ha Une ratios versus [N II]/Ha for each of the galaxies 
in the sample. In the first column, the predicted values of the 
line ratios based on Sokolowski's photoionization model (using 
an absorption-hardened ionizing spectrum and dust-depleted 
abundances in the gas phase) are represented by a series of 
points connected by a thin solid line. The open triangles in 
these panels represent the predictions from Slavin et al.'s dust- 
depleted abundance TML models. Best-fit hybrid models that 
linearly combine the line ratio predictions from pure photoion- 
ization and pure TML models are shown as the thick solid line. 
The best fit is determined visually as the set of points that most 
closely match the data points in all three plots. 

Of the three lefthand panels in Figure 6, the most informa- 
tive is arguably the [O IIIJ/Hq; vs. [N IIJ/Ha diagram. The 
three TML models are clearly separated in this diagnostic dia- 
gram. In the other two panels, the predictions for the log T = 5.3 
and the log T = 5.5 models are too close together to distinguish 
between the two models. Unfortunately, extended [O III] emis- 
sion has not been detected in many of the galaxies, so our anal- 
ysis of the hybrid photoionization/TML models based on the 
[O III]/Ha ratio is limited. The extraplanar |0 III|/Ha profiles 
in the five galaxies in which [O III] has been detected fall into 
two categories: three objects show an increasing [O III]/Ha line 
ratio with increasing |z|, while the [O III]/Hq! ratio drops with 
height in the other two galaxies. The galaxies in the first cate- 
gory tend to be better fit by a TML model with an intermediate 
mixing temperature (log T = 5.3), while those in the second cat- 
egory are better fit by a model with lower mixing temperatures 
(log T = 5.0). 

For five of the galaxies in our sample, we find that the influ- 
ence of the TML regions on the emission-line spectra increases 
with height, such that the contribution to the observed Une ra- 
tios {notto the line flux) increase roughly from ~ 30 to ~ 75%. 
These results are consistent with models in which turbulent 
mixing occurs at higher elevations, where superbubbles break 
out of the thin disk layer, or at the locations where cooling halo 
gas is mixing with the ejected gas. There are exceptions to this 
general rule, however. In UGC 2092 and NGC 4013, the TML 
model contributes to the observed line ratios a near constant 
amount (70% and 40%, respectively), while in NGC 4217, the 
importance of the TML model appears to decrease with height 
(from 40% to 20%). Finally, one of the galaxies in our sample 
(NGC 4302) demonstrates only little deviations from the pho- 
toionization model, and these deviations are not weU explained 
by any of the TML models. 

4.3.2. Shocks 

ShuU and McKee (1979) model interstellar radiating shocks 
such as those due to supemovae events. They include the ion- 
izing effect of the UV precursor in the preshock gas in an effort 



to make their models self-consistent. In their study, Shull and 
McKee discuss ten models. In seven of the models they choose 
as their standard preshock conditions a hydrogen particle den- 
sity of 10 cm~^, a magnetic field Box = 1/xG perpendicular to 
the flow, and cosmic metal abundances. They vary the shock 
velocity from 40 to 130 km s"'. In the other three, they fix 
the shock velocity at 100 km s~' and vary, in turn, the density, 
magnetic field, and metal abundances. 

In the middle column of Figure 6, the predicted line ratios 
from the shock models of Shull & McKee are compared with 
the predictions from the photoionization model (Sokolowski 
1994). The solar abundance shock models are represented by 
the open triangles (the soUd triangle represents Shull and Mc- 
Kee's shock model with depleted abundances), while the pho- 
toionization models are represented by points joined by a thin 
solid fine. The hybrid photoionization/shock model that best 
fits the data is shown as the dark solid line. The models of 
Shull & McKee with shock velocities near 100 km s"' do well 
in predicting the high [O 11I]/Hq; ratios detected in most eDIG 
and the relatively small [O I]/Ha ratios. The biggest problem 
arises when trying to model the galaxies which show a decreas- 
ing [O III]/Ha Une ratio with increasing |z| (recall that two out 
of five galaxies with detected extraplanar [O III] show this be- 
havior). In these galaxies, the [O III]/Hq; ratio drops below that 
which is predicted by the pure photoionization model. Since 
all of the shock models predict enhanced [O Ill]/Ha Une ra- 
tios, they are unable to explain the negative [O lll]/Ha vertical 
gradients. 

In galaxies without extraplanar [O III] data, we have to rely 
heavily on the [S n]/Ha vs. [N Ilj/Ha plot to determine the best 
fitting models. The Une ratios for NGC 3628 and NGC 4013 
can hardly be explained by a combination of shocks and pho- 
toionization (shocks can contribute at most ^ 10% to the ob- 
served line ratios). In the other galaxies, the best-fitting model 
appears to be the depleted abundances model with a shock ve- 
locity of 100 km s"^ In general, the influence of the shocks 
on the observed line ratios appears to increase with increasing 
heights (from ^ 10 - 20% near the disk of the galaxies up to ^ 
40 - 60% at higher |z|). In at least three of the galaxies, the hy- 
brid photoionization/shock model seems to imply that the ion- 
ization parameter decreases from a value of log U = -3 near the 
plane of the galaxy down to log U = -4at heights |z| > 1 kpc. 

4.3.3. Shock H- Precursor 

In a pair of papers, Dopita and Sutherland (1995, 1996) 
present a grid of models of low-density, high-velocity photoion- 
izing radiative shocks. Dopita and Sutherland model photons 
which propagate both upstream through the preshock gas as 
well as downstream into the recombination region of the shock. 
Their grid spans a broad range in shock velocity and magnetic 

parameter, B/nl^^. They assume solar abundances and a low- 
density radiative steady flow shock. The importance of the 
magnetic parameter comes into play because they assume that 
the magnetic field is frozen into the flow. Dopita and Sutherland 
present models both with and without precursor regions. 

In these models, an increase in the shock velocity causes 
an increase in the ionization parameter in the pre- shock gas. 
Since the postshock plasma is the source of ionizing photons 
in radiative shocks, the faster the shock the larger the flux of 
ionizing photons, and hence the higher the ionizing param- 
eter. Dopita and Sutherland also found that an increase in 
the magnetic parameter increases the ionization parameter in 
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the photoionization-recombination zone of the shock, since the 
electron density is inversely proportional to the magnetic pa- 
rameter, and the ionization parameter is inversely proportional 
to the electron density. 

The predictions from Dopita and Sutherland's shock models 
are shown in the right column of Figure 6. The grid spans a 
range in shock velocity from 100 to 500 km s~' (the 500 km s~' 
line is represented by a dashed line) and a range in the magnetic 

parameter B/nl^^ from to 4 fiG cnr'^^ (the 4 /iG cm-*/^ line is 
represented by a dotted Une). Due to the large parameter space 
represented by the Dopita and Sutherland shock models, it is 
difficult to say with much certainty the exact parameters that 
are needed to reproduce the observed line ratios. In general, 
it appears that models combining photoionization and shocks 
with no precursor do better in reproducing the line ratios than 
those using the shocks + precursor models. The measured ra- 
tios also seem to indicate an increase in the shock velocity with 
increasing height. 

5. SUMMARY 

Deep long-slit spectra reaching flux levels of a few 10"'^ erg 
s"' cm"^ arcsec"^ were obtained of nine nearby, edge-on spiral 
galaxies. Ha, [N II] A6583, and [S II] AA6716, 6731 are de- 
tected out to a few kpc in all of these galaxies. Several other 
fainter diagnostic lines are also detected over a slightly smaller 
scale. The relative strengths, centroids, and widths of the var- 
ious emission lines provide constraints on the electron density, 
temperature, reddening, kinematics, and possible source(s) of 
ionization of the eDlG. The main conclusions from the analysis 
are the foUowings: 

• Seven of the nine galaxies in the sample show a general 
increase in the [N II]/Ha and [S II]/Ha line ratios with 
increasing height. The extraplanar [N ll]/Ha line ratios 
reach values in excess of 1.5 in some galaxies. Extra- 
planar [O 111] line emission has been reliably detected 
in five galaxies, and three of those show an increasing 
[O III]/Ha line ratio with height. Extraplanar [O 1] line 
emission has been detected in three of the galaxies, and 
in each case, the [O I]/Ha line ratio increases with in- 
creasing |z|. These trends in the line ratios are similar to 
those observed in other galaxies including our own. 

• He I A5876 is detected in the midplane of four of the 
galaxies, with values ranging from 0.018 to 0.052. This 
suggests an effective stellar temperature range of 36,500 
- 38,500 K and an upper limit to the stellar mass function 
of 42 - 54 M0. Extraplanar He I emission is detected in 
only one galaxy, NGC 2820. The He 1/Ha and [N ll]/Ha 
line ratios in the eDIG of this galaxy are consistent with 
O-star photoionization, although the other ratios in this 
galaxy suggest the presence of an additional source of 
ionization. 

• For all but one galaxy in the sample (NGC 4302), pho- 
toionization by massive OB stars does not appear to be 
sufficient to explain the fine ratios in the eDIG. At least 
one other source of ionization appears to be needed. Hy- 
brid models that combine photoionization by OB stars 
and photoionization by turbulent mixing layers (TMLs) 
or shocks provide the best fit to the data. In contrast 
to the pure photoionization models, these hybrid mod- 
els are able to reproduce the run in the [O III]/Hq! line 



ratios observed in many of the galaxies. Overall, the 
hybrid photoionization/TML models do a better job of 
explaining the observed line ratios than the photoion- 
ization/shock models, although in many cases the pho- 
toionization/shock models do almost or just as well. The 
contribution of the turbulent mixing layers (or shocks) to 
the observed line ratios appears to increase with increas- 
ing height. These results are consistent with models in 
which turbulent mixing and shocks occur at higher ele- 
vations, where superbubbles break out of the thin disk 
layer, or at the locations where cooling halo gas is mix- 
ing with the ejected gas. 

• Three galaxies in the sample appear to show significant 
vertical velocity gradients. In NGC 3628, the presence 
of a prominent dust lane is affecting the velocity mea- 
surements and producing large variations in excess of 
100 km s~^ within a range in height of only 0.5 kpc. 
Monotonic velocity gradients of order 50 - 70 km s~' 
kpc"' appear to exist in NGC 2820 and NGC 4013 (a 
gradient may also be present in NGC 4302, but the un- 
certainties are large). These gradients can be explained 
if the rotational velocities of the eDIG are lower in the 
eDIG than in the disk, as predicted by the galactic foun- 
tain model. The upper limits on the velocity gradients in 
the other galaxies (~ 30 - 50 km s~' depending on the 
galaxy) are not inconsistent with this model. 
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Appendix: Notes on Individual Objects 



UGC 2092 

This is a highly inclined (i = 86°; Guthrie 1992) Scd galaxy 
that lies at a distance of 72 Mpc (PBS). PBS detected a bubble- 
like structure approximately 20" NE in the plane of the galaxy. 
We positioned the slit so that it passes directly through this 
structure (Fig. 1). UGC 2092 is the most distant galaxy in the 
sample. In order to detect even the strongest emission lines, the 
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data were spatially binned by a factor of two, yielding a spatial 
scale of ~ 1.1 kpc per spatial element. The TML models have 
difficulties explaining the observed line ratios. Higher mixing 
gas temperatures seem to be favored, but if this were the case, 
then the TML model would have to account for almost all of 
the observed emission on the west side of the galaxy, while the 
east side would require an ionization parameter smaller than 
-5.0. The shock models do a better job of explaining the ob- 
served line ratios. A shock velocity of 100 km sec~^ and de- 
pleted abundances appear to best model the shock environment. 
According to the [S II]/Ha vs. [N II]/Ha plot, shocks play a 
more important role on the east side of the galaxy, contribut- 
ing about 60% to the observed line ratios, while on the west 
side, only 15% - 45% of the observed Une ratios is produced 
through shocks. Dopita & Sutherland's shock models also sug- 
gest a stronger influence from shocks on the eastern side of the 
galaxy. The observed line ratios can best be modeled by emis- 
sion originating in the recombination region of the shock, with 
a magnetic parameter of ~ /uG cm^/-^. 



UGC 3326 

UGC 3326 is a Scd galaxy which lies at a distance of 48 Mpc 
(PBS) with an incUnation angle of 90° (Guthrie 1992). PBS 
detected a faint plume on the SW side of the galaxy. The plume 
extends outward at an angle of 52° away from the axis of the 
galaxy out to a height ^ 1.3 kpc. At the end of the plume is 
a relatively bright cloudlike structure. The slit was centered 
on the galaxy and ran directly through this plume (Fig. 1). 
UGC 3326 also required 2x spatial binning in order to detect 
even the strongest emission lines, yielding a spatial scale of 
~ 700 pc per spatial element. As in the case of UGC 2092, 
only [N II]/Ha and [S II]/Ha are available to estimate the in- 
fluence of a secondary source of ionization. The observed line 
ratios near the disk of the galaxy (out to ~ 700 pc) are consis- 
tent with the photoionization model and do not require a sec- 
ondary ionization source. At higher |z| (on the northwest side 
of the galaxy), however, the [S II]/Ho; line ratios are larger than 
what can solely be explained by the photoionization model. The 
TML model suggests that at higher elevations, about 40 - 80% 
of the line ratios is produced by the mixing layer gas, while 
the shock model suggests that 20 - 40% of the line ratio arises 
from shocks (based on the 100 km sec~\ depleted abundances 
model). Unfortunately, with only Ha, [N 11] and [S 11] lines 
observed in this galaxy, we cannot further constrain the data or 
verify which model best explains the observed eDIG emission. 



UGC 4278 

R96 imaged UGC 4278, an SB(s)d galaxy with an inclina- 
tion of 90° (TuUy 1988), and detected prominent plumes on 
both sides of the nucleus of the galaxy, out to a height of ^ 
1.4 kpc. R96 suggests that this may trace outflow from a nu- 
clear starburst. Our sUt was positioned so it passes through the 
nucleus of the galaxy, and therefore through the emission-line 
regions detected by R96 (Fig. 1). UGC 4278 lies at a distance 
of 10.6 Mpc (Tully 1988), yielding a spatial scale of ^ 160 pc 
pixeL^ The [O III]/Ha data in UGC 4278 tend to compUcate 
the identification of the secondary source of ionization. While 



the data for this line ratio, as well as [S II]/Ha, suggest the need 
for an additional source of ionization, they conflict in the deter- 
mination of exactly what that secondary source might be. While 
the [S II]/Ha line ratio shows a general increase with height, 
the [O 111]/Hq; line ratio show a general decrease with height. 
While the [S II]/Ha vs. [N n]/HQ; data points may be explained 
with a photoionization model with ionization parameters from 
~ -2.0 to -3.0 mixed with shocked gas with a shock velocity 
of 100 km sec"^ in a depleted abundance region, neither this 
model nor any other shock model can explain the decrease in 
the [O lll]/Ha line ratio with increasing height. Similar prob- 
lems occur with the TML models. While the TML models with 
low mixing gas temperatures (such as log T = 5.0) have neg- 
ligible [O 111] emission, and so could account for a decrease 
in this line relative to Ha if TML becomes more important at 
large heights, the predicted [S II]/Ha and [N ll]/Ha line ra- 
tios for this model lie close to the photoionization curve, and 
so limits the available parameter space. Looking at Figure 6, 
the observed [S Il]/Ha and [N ll]/Ha ratios are clearly larger 
than predicted by any combinations of O-star photoionization 
and log T = 5.0 TML models. 

The remarkably small [N II]/Ha line ratios in the disk of 
this galaxy have been noted previously by Goad & Roberts 
(1981). In their study of superthin galaxies, one of which was 
UGC 4278, fliey derive a mean [N IIJ/Ho; line ratio of 0.12 for 
this object, consistent with our measurements. They conclude 
that the small | N 11 |/Ha ratios are a result of nitrogen deficiency 
in this galaxy. This suggests that star formation proceeds on a 
much longer time scale than in most galaxies, since nitrogen 
is a secondary product of nucleosynthesis. It also implies that 
large-scale shocks are rare in UGC 4278. 



NGC 2820 

NGC 2820 lies at a distance of 20 Mpc, yielding a spatial 
scale of ^ 300 pc pixel"' . The position of the slit was selected 
based on the narrowband images reported in Paper 1. As dis- 
cussed in §3.1.2, NGC 2820 is the only galaxy in our sample 
where extraplanar He 1/Ha emission has been detected out to 
|z| ~ 1 kpc. It is also the only galaxy in which all four main 
line ratios are detected out to large vertical heights, therefore 
providing strong constraints on the nature of secondary ioniza- 
tion source. The main concern with NGC 2820 is that the ma- 
jority of the [O lll]/Ha vs. [N 11]/Hq data points lie below the 
photoionization curve. The TML models with low mixing gas 
temperatures (log T = 5.0) also have low |0 lll]/Ha values. Us- 
ing this model in conjunction with the photoionization model, 
the [O lll]/Ha and [O 1]/Ha data suggest that turbulent mixing 
layers contribute roughly 35 - 55% to the observed line ratios. 
Unfortunately, the [S II]/Ha data do not support this conclu- 
sion. The [S n]/Ha data suggest a higher mixing gas tempera- 
ture (log T = 5.3), with the contribution to the line ratios aris- 
ing from this region apparently increasing with heights from 
20 to 80%. The comparison with the photoionization/shock 
hybrid model suggests that the data are best fit by the model 
with shock velocity = 100 km sec"^ and depleted abundance. 
With the exception of those [O lll]/Ha data points that fall be- 
low the photoionization curve, the remaining data suggest that 
the ionization parameter decreases from ~ -3.0 to -4.0 while 
the contribution of shocks to the line ratios rises from ~ 10 to 
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40%. The Dopita & Sutherland models are perhaps more suc- 
cessful at explaining the observed trends in the line ratios. The 
[N II]/Ha, [S II]/Ha, and [O III]/Ha Une ratios lie in the region 
of the (shock + precursor) model with a magnetic parameter of 
~ /iG cm3/2. 



NGC 3628 

NGC 3628 is a starbursting Sb galaxy which is a member 
of the Leo Triplet. Located at a distance of 6.7 Mpc, it has 
an inclination angle of 87° (Tully 1988). Imaged by Fabbiano, 
Heckman, and & Keel (1990), they detected a large plume ex- 
tending about 9 kpc from the west side of the galaxy (see also 
Dahlem et al. 1996). It is not exactly perpendicular to the 
disk of the galaxy, but lies at a position angle of 210° (while 
the disk of the galaxy lies at a position angle of 104°). In or- 
der to cover as much of this plume as possible, we positioned 
the slit to extend only to vertical heights to the south side of 
the disk, following the extent of the plume, rather than cen- 
tered on the disk of the galaxy (see Fig. 1). Fabbiano et al. 
also obtained spectroscopic information along the minor axis 
of NGC 3628, though not through the detected plume. They 
detected a midplane [N n]/Ha ratio of 0.4 which increased to 
> 1.0 at higher |z|. [S II]/Ha was also detected and observed to 
increase with |z|, having a midplane value of 0.3 and rising to 
^ 0.8 at large heights. Shock ionization was proposed by Fab- 
biano et al. (1990) to explain these line ratios. Figure 6 shows 
that pure shock or photoionization/shock hybrid models are un- 
able to reproduce the small [O in]/HQ line ratios in the extra- 
planar plume; this is different from the conclusion of Fabbiano 
et al (1990) who did not have constraints on the [O III] emis- 
sion. The photoionization/TML hybrid models with a mixing 
gas temperature of log T = 5.0 are more successful at explaining 
the observations. Based on this model, the TML region would 
contribute roughly 50% to the observed line ratios, although the 
exact amount differs between the [O lll]/Ha and [S II]/Hq! data 
points. The plume may represent the interface where hot wind 
gas mixes with entrained disk material (Dahlem et al. 1996). 



NGC 4013 

NGC 4013 is a Sb galaxy at a distance of 17 Mpc (260 pc 
pixel"'; Tully 1988) and with an inclination of 90° (Bottema 
1995). R96 detects extraplanar gas features extending up to |z| 
^ 2.5 kpc on the northeast side, as well as a less prominent 
line-emitting region on the southwest side. Our recent imaging 
study (Paper 1) confirms the presence of eDIG in this object. 
We positioned our slit to cover the northeast region. The pho- 
toioni/ation/TML hybrid model with log T = 5.0 works well for 
NGC 4013. While there are only a few [O lll]/Ha data points, 
they are consistent with the |S ll]/Ha data points, suggesting 
an ionization parameter of ^ -3.8 and implying that approxi- 
mately 50% of the Une ratios originate from the TML region. 
The photoionization/shock hybrid model is unable to reproduce 
the observed emission-line ratios. 



NGC 4217 



NGC 4217 is a Sb galaxy at a distance of 17 Mpc (260 pc 
pixel"'; Tully 1988) with an inclination of 86° (R96). R96 de- 
tects two bright plumes of extraplanar emission extending ~ 
2 kpc from either side of the disk on the southwest side of 
NGC 4217. Our slit was positioned so that it passes through 
both plumes of extraplanar gas. [O 1] was detected out to ~ 
1.3 kpc on either side of NGC 4217. Based on the [O I]/Ha 
line ratios, in addition to the observed [S II]/Ha line ratios, we 
find that the photoionization/TML hybrid models predict that 
the ionization parameter decreases from -4.0 to -5.0, while the 
influence of TMLs on the observed Une ratios remains relatively 
constant at ~ 40%. The best fit photoionization/shock model 
suggests an ionization parameter which decreases from -3.5 to 
^.8, and a shock model with a velocity of 100 km sec"' and 
depleted abundances, which provides only 10 to 20% of the 
observed line ratios. According to the Dopita & Sutherland 
shock models, the observed emission is inconsistent with shock 
+ precursor emission. Although there is significant scatter in 
the data, the observed line ratios appear to be well described 
by emission originating from a region with a magnetic param- 
eter < 1 /xG cm^/^ and a shock velocity which increases with 
increasing height (up to almost 400 km sec"'). 



NGC 4302 

NGC 4302 was imaged by both R96 and PBS. It is an Sc 
galaxy at a distance of 16.8 Mpc (250 pc pixel"'; Tully 1988) 
with an inclination of 90° (R96). PBS detected an extremely 
faint emission-Une structure emerging 0.73 kpc from the plane 
close to the nucleus of the galaxy. R96 detects faint, diffuse 
extraplanar emission at galactocentric radii < 4 kpc up to |z| ^ 
2 kpc. Therefore, we chose to place our sUt through the cen- 
ter of the galaxy to incorporate both structures. Unfortunately, 
the uncertainties on many of the extraplanar line ratios are large 
and prevent us from making strong statements on the source of 
ionization in the eDIG of NGC 4302. All of the data points 
are consistent with the photoionization model, within the large 
measurement errors. 

NGC 4302 was also observed spectroscopically by Collins & 
Rand (2001). The position of their slit is 20" southeast of ours, 
or given a distance of 16.8 Mpc, ^1.6 kpc southeast. They de- 
tect the prominent [N 11] and [S II] lines up to z = 2 kpc on either 
side of the galaxy, but poor weather conditions prohibited the 
detection of [O 111]. Their [N II]/Hq! line ratios increase from 
0.4 up to 1.4 at |z| = 1.5 kpc, and their detected [S 1I]/[N 11] 
Une ratios remain constant at 0.6. These values suggest a gas 
temperature ranging from 6600 to 10,000 K. While their line 
ratios differ from ours, we do see the same general increase in 
the emission Une ratios with Izl . 



NGC 5777 

In NGC 5777, an Sb galaxy with an inclination of 83° 
(Guthrie 1992) located at a distance of 25 Mpc (380 pc pixel"'; 
PBS), PBS detected a prominent extraplanar plume extending 
1.1 kpc directly north-east from a bright H II region located on 
the southeast side of the galaxy. We positioned our slit so that it 
was centered on the galaxy disk and passed through this plume. 
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A photoionization/TML hybrid model with a mixing gas tem- 
perature of log r = 5.3 best fits the data. This fit suggests that 
the ionization parameter decreases from ^.0 to -5.0 while the 
contribution from the TML region to the emission Une ratios 
increases from ^ 20 to 80%. The photoionization/shock hybrid 
model is less successful at explaining the data. The best fit is 
found with a model with a shock velocity of 100 km sec"' and 
depleted abundances, where the contribution from the shock 



models to the emission line ratios increases from ^ 15 to 40% 
with increasing heights. The Dopita & Sutherland models also 
have difficulties reproducing the data. The [O III]/Ha ratios 
measured at high |z| favor a shock + precursor model, while the 
[S 11]/Hq! data suggest that the emission originates from shocks 
without precursor. These ratios also suggest different values for 
the magnetic parameter. 
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Fig. 1 .— Position of tlie slit superposed on tlie continuum-subtracted Hq image of (a) tlie NE side of UGC 2092, (b) tlie SW side of UGC 3326, (c) UGC 4278, 
(d) NGC 2820, (e) NGC 3628, (f) NGC 4013, (g) NGC 4217, (h) NGC 4302 and (i) NGC 5777 (see Table 1 for references to the images). The position of the sUt is 
chosen to pass through at least one known region of eDIG. In most cases the slit is positioned perpendicular to the plane of the galaxy. In the others, it is positioned 
to maximize coverage of the eDIG. UGC 2092 is rotated clockwise by 32° such that the disk of the galaxy runs vertically, while UGC 3326 was rotated clockwise 
by 84° . The other figures have not been rotated (North is at the top of the figure and east to the left). 

Fig. 2. — Emission line spectra as a function of position along the slit. For display purposes, the spectra are scaled individually (or in sets of two or three) with 
the scaling factor in relation to the spectrum containing the Ha + [N II] complex listed in the upper right comer of each spectrum. The expected redshifted positions 
of the emission Unes as determined from the systemic velocity of the galaxies are indicated by the dashed lines. Tickmarks on the horizontal axis are separated by 
25 A in the observer's restframe. In the cases where no line is detected, a representative spectrum from near the disk of the galaxy is presented. The presence of the 
Na ID AA5897.6, 5891.6 absorption Unes near He I is apparent in some of the galaxies (e.g., NGC 3628, NGC 4217, NGC 4302 and UGC 3326). 

Fig. 3.— Vertical profiles of [N H] A6583/Ha, [S H] A6716/Ha, [S H] A6716/[N 11] A6583, [O m] A5007/Ha, [O I] A6300/Ha, He I A5876/Ha, Ha/H/3, [S H] 
A6716/[S n] A6731, and [O HI] A5007/[O HI] A4959. This last ratio is expected to be constant (~ 3) within the uncertainties. The error bars represent one-cr 
uncertainties. Also shown in each panel are the spatial profiles of the Ha emission (solid line) and continuum emission (dashed line), both of which have been 
arbitrarily scaled for display purposes. Only ratios involving detected lines have been included in this figure. 

Fig. 4. — Temperature profile and ionization fraction in the eDIG. In the bottom left panel is a plot of the temperature profile (in units of IC* K) based on the [N II] 
A6583/Ha line ratio presented in the upper left panel [see eqn. (2)]. Using this temperature profile, [S 11] A6716/Ha, [S II] A6716/[N II] A6583, [O 111] A5007/Ha 
and [O I] A6300/Ha have been calculated and superimposed on the observed line ratios [see §4.2 and eqns. (3) through (6) for details]. The lines represent different 
ionization fractions ranging from 12.5% (first solid line) to 100% (last dotted line). The top two panels to the right of the [N II]/Ha plot represent the S+/S ratio; the 
bottom, middle plot represents the O^/O ratio; and the bottom, right plot represents the H+/H ratio. The [N II]/Ha ratios cannot be used to derive the temperature 
profile in objects with strong eDIG [O III]/Ha ratios (> 0.25); four objects are therefore excluded from the analysis (UGC 4278, NGC 2820, NGC 4302, and 
NGC 5777; see text). ^ 

Fig. 5. — Velocity and line width profiles. In three of the panels, the velocity centroids of the Ha, [N II] A6583 and [S II] A6716 emission lines are plotted as a 
function of position along the slit. The velocities are plotted relative to the systemic velocity in km s"' . The dot-dashed line represents the average position of the 
line centroid. One-u error bars have been computed based on the accuracy in detecting the true line centroid as a function of signal-to-noise ratio. The bottom, right 
figure presents the measured Ha line widths in km s~' as a function of vertical height. These measurements have been corrected for the width of the instrumental 
profile. 

Fig. 6.^ (a) Comparisons of observed [S II] A6716/Ha, [O III] A5007/Ha, [O I] A6300/Ha, and [N II] A6583/Ha line ratios with the predictions of pho- 
toionization, turbulent mixing layer, and shock models. The data for UGC 2092 are the filled circles with error bars. The series of points joined by a thin solid 
line in the panels in the left and middle colunms represent the predictions from photoionization models of a hardened spectrum and dust depleted abundances with 
ionization parameters (log U) of -3.0, -3.5, -4.0, -A.2, -4.5, and -5.0 (note that [N ITj/Ha decreases with increasing log U\ Sokolowski 1994). The left column 
overlays the predictions from the turbulent mixing layer (TML) models of Slavin et al. (1993). These values are derived from models with a hot gas velocity of 
25 km s"', depleted abimdances, and an intermediate temperature after mixing of log T = 5.0, log T = 5.3, or log T = 5.5 (represented by open triangles with size 
increasing with temperature). The hybrid photoionization/TML model which best fits the data is shown as a thick soHd Une. This model incorporates TMLs with 
a mixing temperature of log T = 5.5 (see text for more detail). The middle colunm overlays shock model predictions from ShuU & McKee (1979). The values 
plotted are calculated from solar abundance models with shock velocities of 90, 100, and 1 10 km s~^ (represented by open triangles of increasing size). A model 
with a shock velocity of 100 km s~' and depleted abundances (filled triangle) is also shown for comparison. For this galaxy as well as the others, the best fit to the 
photoionization/shock hybrid model is derived using the shock model with depleted abundances. The panels in the right colunm presents the grid of predictions 
from the shock and shock+precursor solar abundance models of Dopita & Sutherland (1995). The shock velocity ranges from 100 to 500 km s~' (the 500 km s~' 

line is represented by a dashed line) and the magnetic parameter B/nl^^ ranges from to 4 fj,G cm^/^ (the 4 fj,G cm^/^ line is represented by a dotted hne). (b) Same 
as Figure 6a but for UGC 3326. The best fit for the photoionization/TML hybrid model uses a mixing temperature of log T = 5.3. (c) Same as Figure 6a but for 
UGC 4278. The best fit for the photoionization/TML model uses a mixing temperature of log T = 5.0. (d) Same as Figure 6a but for NGC 2820. The best fit for the 
photoionization/TML hybrid model uses a mixing temperature of log T = 5.3. (e) Same as Figure 6a but for NGC 3628. The best fit for the photoionization/TML 
hybrid model uses a mixing temperature of log T = 5.0. (f) Same as Figure 6a but for NGC 4013. The best fit for the photoionization/TML hybrid model uses a 
mixing temperature of log T = 5.0. (g) Same as Figure 6a but for NGC 4217. The best fit for the photoionization/TML hybrid model uses a mixing temperature of 
log T = 5.0. (h) Same as Figure 6a but for NGC 4302. These data do not require a secondary source of ionization, (i) Same as Figure 6a but for NGC 5777. The 
best fit for the photoionization/TML hybrid model uses a mixing temperature of log T = 5.3 
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Fig. 1. — See jpg image. 

Fig. 1. — (Cont'd) See jpg image. 

Fig. 1. — (Cont'd) See jpg image. 







[0 III] 


He 


UGC 


2092 

[0 1] 


Ha,[N II] 


[S II] 




1 1 1 1 


1 III III 
1 1 








1 II 1 1 1 
1 


II II 1 


Z Height (kpc) 

Last 

-2-10 1 


1 
1 

1 
1 

1 

1 1 ' 1 


1 1 

1 ' 1 1 ' 








j 1 1 .00 
Ijl ,.00 
Ij 1 1 .00 


ill -I 

Wii'iiiVi,'' 

4^fl 



wavelength (Angstroms) 

Fig. 2.— 



13 



UGC 3326 

H,!? [0 III] He I [0 I] Ha.[N II] [S II] 



CO 












1 1 II 1 1 1 
1 1 1 


III 1 

1 1 


o 

Cl 

^ o 
ai 
a) 

X 
M 


ill 


%i 


fi 


1 
1 

1 


1 
1 

f'^ ,fft?ihl^WV'>''^ 

1 




' ( ' 1 .00 

( ! 1.00 

I'l ! 
' 111 

1 ' !\ 

,ifcy'|;J'l,V.-*vv 


j,'0-58 


1 

^ 1 


1 
1 

1 




1 
1 

1 




1 
1 

1 






i I, 1-=° 
i| 1 1 -00 


1 1 
1 1 



wavelength (Angstroms) 



[0 III] 



UGC 42^8 
He I [0 I] 




- I* 



0.59 



fO 54 



I 1 -02 




\k^i'lV''VjyiV%v^ri 

I ' i' ' I 
1 J I'' ^'■^'i' 



0.04 

0.06 

m>Ajl<,..'iV 



•Jl.'vl 
^ ii; III 



Ha.[N II] [S II] 



I I 1 .00 I I 0..3 
I I l)| III 




u.uz 
0.05 



v;*i I'll 

1 .00 




I ' 1.00 



111 



(l I I ^ 



ll|f(lll"y! 
,1,1 IPl5.5 



0,14 



Si 



.J: 



)| I 0, 1 9 



0,37 

iii''iM' 

I 1 0,48 



[0 III] 



NGC 2820 
He I [0 I] 



il 0.26 



) e.29 



! 0.3; 



0.3^^ 



I'l C'.2' 




,1 0.43 



I 0.38 



Pl/l^iW 



0.07 

4>Vi,iii,*,V 

0.06 



0,28 



Ha,[N II] [S 

,1^! iHiWiJ'-*j ,>f (''1 i'|iri«''( I 



0,07 



0,05 



j 0,03 
0,02 



i,,ii;vf'j'Vf% 

0,04 




I I 



1,00 I ;o,33 

1 ,00 



.vvviv^luS 

J 1 1 ,00 



( )0,1 1 

III )l^ 

1 10,18 



wavelength (Angstroms) 

Fig. 2.— (Cont'd.) 



Wavelength (Angstroms) 



14 



[0 III] 



NGC 3628 

He I [0 I] 



Ha,[N II] [S II] 



'l 0.25 



j 0.28' 




H/3 



[0 III] 



NGC 4013 
He I [0 I] 





[ 0.24 



'fir 



1 .01 
0.24 
0.07 
0.10 



■ ■ li.oo 

1.00 



Ha,[N II] [S 



wavelength (Angstroms) 

Fig. 2.— (Cont'd.) 



wavelength (Angstroms) 



15 




NGC 5777 NGC 5777 




wavelength (Angstroms) 

Fig. 2.— (Cont'd.) 



wavelength (Angstroms) 



16 




17 




18 




Fig. 3.— (Cont'd.) 



20 




Fig. 3.— (Cont'd.) 



21 



n 




SH I S'O Sl'O l-'O SO'O 8 9 * Z 

oRoy sun 

Fig. 3.— (Cont'd.) 



22 




23 




24 




OBoy sun 

Fig. 3.— (Cont'd.) 



25 



Table 1 
Sample 



Galaxy 


R.A. (12000)" 


Dec (J2000)" 


0,5" 


Disk P.A." 


Morphological 


Dist." 


Incl." 


Reference for 




hh mm ss 


dd mm ss 


(') 


(°) 


Type" 


(Mpc) 


(°) 


Imaging Data 


UGC 2092 


02 36 30.0 


+07 18 00 


3.16 


32 


Scd 


72' 


86" 


Pildisetal. 1994 


UGC 3326 


05 39 36.0 


+77 18 00 


3.55 


84 


Scd 


48' 


90" 


Pildis et al. 1994 


UGC 4278 


08 13 59.0 


+45 44 43 


4.68 


172 


SB(s)d 


10.6 


90 


Rand 1996 


NGC 2820 


09 21 47.1 


+64 15 29 


4.07 


59 


SB(s)c 


20.0' 


90 


Paper I 


NGC 3628 


1120 16.3 


+13 35 22 


14.79 


104 


Sb 


6.7 


87 


Fabbiano et al. 1990 


NGC 4013 


1158 31.7 


+43 56 48 


5.25 


66 


Sb 


17 


90'' 


Rand 1996, Paper I 


NGC 4217 


12 15 50.9 


+47 05 32 


5.25 


50 


Sb 


17 


86" 


Rand 1996 


NGC 4302 


12 21 42.5 


+ 14 36 05 


5.50 


178 


Sc 


16.8 


90" 


Rand 1996 


NGC 5777 


14 51 18.3 


+58 58 35 


3.09 


144 


Sb 


25' 


83" 


Pildis et al. 1994 



"Values taken from de Vaucouleurs et al. 1991. 
''Values taken from Tully 1988 unless otherwise noted. 

•^References for distances: UGC 2092, UGC 3326, & NGC 5777 (Pildis et al. 1994), NGC 2820 (Hummel & van der Hulst 1989). 

''References for inclinations: UGC 2092 & UGC 3326 (Guthrie 1992), NGC 4013 (Bottema 1995), NGC 4217 & NGC 4302 (Rand 1996), NGC 5777 (Guthrie 1992). 



Table 2 
Observing Logs 



Galaxy 


Total Exposure 


Number of 


Slit 


Distance from 




Time 


Exposures 


RA. (°) 


Nucleus 


UGC 2092 


5.75 hours 


9 


122 


20"NE 


UGC 3326 


4.5 hours 


7 


107 


16"SW 


UGC 4278 


5.0 hours 


8 


83 


0" 


NGC 2820 


5.0 hours 


7 


139 


19"NE 


NGC 3628 


3.75 hours 


5 


210 


20"W 


NGC 4013 


4.5 hours 


6 


155 


40"NE 


NGC 4217 


4.5 hours 


6 


144 


50"SW 


NGC 4302 


4.5 hours 


6 


88 


0" 


NGC 5777 


4.5 hours 


6 


54 


50"SE 



Table 3 

Reddening Factors for Important Line Ratios if Av = 1 



Line Ratio 


Reddening 


[NH] A6583/Ha 


1.004 


[S II] A6716/Ha 


1.022 


[O I] A6300/HQ; 


0.963 


He 1 A5876/Ha 


0.898 


[O III] A5007/Ha 


0.751 


Ha/H/3 


1.381 


[S II] A6716/[NII] A6583 


1.021 


[S II] A6716/[S II] A6731 


0.998 


[O ni] A5007/[O m] A4959 


1.012 
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